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The creatine kinase (CK) reaction plays a critical role in skeletal muscle function, and can be studied 
non-invasively using phosphorus (^^P) saturation transfer (ST) techniques. However, due to the low MR 
sensitivity of the ^^P nucleus, most studies on clinically approved magnetic fields (<3.0 T) have been 
performed with coarse resolution and limited tissue coverage. However, such methods are not able to detect 
spatially resolved metabolic heterogeneities, which may be important in diseases of the skeletal muscle. In 
this study, our aim was to develop and implement a ^^P-MRI method for mapping the kinetics of the CK 
reaction, and the unidirectional phosphocreatine (PCr) to adenosine triphosphate (ATP) metabolic fluxes 
in muscles of the lower leg on a clinical 3.0 T MR scanner. We imaged the lower leg muscles of ten healthy 
volunteers (total experimental time: 40 min, nominal voxel sizes 0.5 mL), and found statistically significant 
differences between the kinetics of the CK reaction among muscle groups. Our developed technique may 
allow in the future the early detection of focal metabolic abnormalities in diseases that affect the function of 
the skeletal muscle. 

Phosphorus (^^P) magnetic resonance (MR) is a unique non-invasive tool for studying muscle physiology, 
and has been used extensively for studying diseases that adversely affect muscle function, such as heart 
failure\ stroke, congenital myopathies^, and ischemic heart disease^. A unique feature of saturation transfer 
(ST) ^^P-MR techniques is their ability to measure the kinetics of important metabolic reactions that involve 
phosphorus containing metabolites^. High- energy phosphate is reversibly transferred between adenosine tripho- 
sphate (ATP) and phosphocreatine (PCr) through the creatine kinase (CK) reaction. Disturbances in the kinetics 
of the CK reaction are found in many disease that affect myocardial energy metabolism^'^. However, in the skeletal 
muscle the physiological significance of alterations on the CK exchange rate are less understood. 

The kinetics of the CK reaction have been studied extensively using surface radiofrequency (RF) coils and, 
either unlocalized^ or localized methods with limited tissue coverage^. However, such methods lack the ability to 
measure muscle- specific energy metabolism and therefore cannot identif)^ potential heterogeneneities in muscle 
function in health and disease. The development of relatively high-resolution and large tissue coverage imaging 
approaches for identifying spatial heterogeneities in the kinetics of the CK reaction may enable a better under- 
standing of the CK role in the muscle. 

A major roadblock for the development of high-resolution ^^P-MR methods with sufficient tissue coverage is 
the low MR sensitivity of the ^^P nucleus and the relatively low concentration of ^^P containing metabolites in the 
human tissue. The increased availability of ultra-high field (UHF) whole body magnets (>7.0 T) has allowed the 
development of time efficient ^^P-MR methods that provide adequate tissue coverage and significantly increased 
spatial resolution^"^\ However, only a few tens of UHF systems are currently operational worldwide^^, while most 
of the clinical MRI work is performed almost exclusively in fields of 3.0 T or lower. Therefore, it is equally 
important to develop ^^P-MR tools for the diagnosis of complications in the skeletal muscle in fields of 3.0 T or 
lower. 

In this study, we aimed at developing and implementing a three-dimensional saturation transfer (ST) ^^P-MRI 
method capable of mapping the CK reactions kinetics in the entire muscle of the lower leg on a clinically approved 
3.0 T whole body magnet. Our method maps the unidirectional conversion rate of PCr-to-ATP at relatively high 
spatial resolution (i.e. 0.5 mL voxel size), in a total acquisition time of 40 min, and can potentially be used in a 
routine clinical setting. 

Results 

In order to estimate the pseudo first-order forward rate constant (kf) of the CK reaction, we need to measure the 
phosphocreatine (PCr) signal while we saturate the y-adenosine triphosphate (y-ATP) resonance for different 
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Figure 1 | Unlocalized ST-^^P-MRS of the lower leg. A) Fully relaxed ^^P spectrum without ST preparation (^at 
complete saturation (arrow) of the y-ATP resonance (tsat = 6.84 s) (right). B) Series of ST spectra at different tsa 
(Pearson's product moment correlation = 0.9992). 



= 0 s) (left) and '^P-ST Spectrum with 
, C) PCr signal intensity fitted to Eq.2, 



durations (i.e. the progressive ST experiment) Assuming a two- 
pool exchange system between PCr and y-ATP and complete sat- 
uration of y-ATP, we can estimate the exchange rate between the two 
metabolites by solving the modified, for chemical exchange, Bloch 
equation^^ (see methods). To confirm the efficiency of our saturation 
pulses, we acquired unlocalized ^^P spectra in the entire volume of the 
lower leg muscles in all of our volunteers. Typical ^^P spectra can be 
seen in Fig. 1, where the full width at half maximum of the PCr 
resonance peak was 12.3 ± 4.8 Hz (mean ± SD) across all subjects. 
As we show in the same figure, the ST module saturates y-ATP to 
levels not measurable above the noise. 

Fat infiltration, especially in high body mass index (BMI) subjects, 
can reduce the volume fraction of lean muscle in a certain volume of 
tissue affecting quantification of PCr concentration. In order to 
account for the fat content in the muscle, the Hierarchical IDEAL 
method^^ was applied across all subjects. The fat volume fractions in 
the TA were 0.035 ± 0.006 across all subjects, significantly lower 
than in the GL (0.056 ± 0.019, P = 0.012), the GM (0.051 ± 0.020, P 
< 0.001) and the S (0.063 ± 0.020, P = 0.002) muscles. Figure 2 
shows a case of a higher BMI (30.4) subject. The water/fat separation 
method shows fat infiltration in the leg muscles. The fat volume 
fraction in the GM was 0.11, and 0.04 in the TA, which, if not 
accounted for (especially the GM), can lead to under- estimation of 
the PCr concentration. 

In the turbo spin echo (TSE) imaging method, the echo amplitude 
is modulated as a function of the echo position in /c-space, resulting 
in blurring^^. We estimated the point spread function (PSF) of our 
imaging method as the result of the echo modulation to be at 1.6 
pixels in the anterior to posterior direction, where the echo train is 
sampled. 

Figure 3A shows an anatomical cross-section of the lower leg 
muscles of a lean subject (BMI = 19.5) and a cross-section of the 
PCr concentration map (Fig. 3B). PCr in the same cross-section 
decreases with increasing saturation time (tsat) as shown in Fig. 3C. 
By segmenting signals in the PCr images and fitting data to Eq.2 



(methods), we estimated /cf and the unidirectional flux of PCr to form 
ATP (Vf) for four muscle groups of the leg [Gastrocnemius Lateralis 
(GL) and Mediahs (GM), Soleus (S), and Tibiahs Anterior (TA)]. We 
estimated the intrinsic spin lattice relaxation (Ti ) in the presence of 
saturating irradiation^^ from Eq.3. The results are summarized in 
Table 1. In the TA, kf was significantly lower than the GL (P = 
0.002), the GM (P = 0.029) and the S (P = 0.037). The metabolic 
fluxes, Vf, in the TA were lower than both the GL (P < 0.001) and the 
GM (P = 0.043). We did not find any statistically significant differ- 
ences among values in the four muscle groups. 

Discussion 

Our study focused on the development and implementation of a ST - 
^^P-MRI method capable of obtaining muscle specific measurements 
of the CK reaction and metabolic fluxes of PCr to ATP in the entire 
volume of the muscles of the lower leg using a clinical high-field 
(3.0 T) scanner. Our method enables such measurements with rela- 
tively high spatial resolution (i.e. 0.5 mL voxel size) and large muscle 
coverage within an acquisition time of 40 min. 

Our imaging approach (Fig. 4B) is based on a fast 3D-TSE sequence, 
which excites and acquires a single resonance of the ^^P-MR spectrum 
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Figure 2 | Water/fat fractions in the lower leg. The water (left) and fat 
(right) images of a subject (BMI = 30.4) are shown. The fat fraction (right) 
shows fat infiltration in the leg muscles. 
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Figure 3 | Saturation transfer ^^P-MRI. A) Anatomical proton image 
where four major muscle groups are identified. B) PCr concentration map 
in a cross-section of the lower leg. C) ST-^^P-MRI of the same cross-section 
at different saturation times. PCr signal intensity decreases as the system is 
allowed to exchange for longer time. 



(i.e. PCr). By trading spectral information for imaging speed, we are 
able to accelerate acquisition time and increase spatial resolution 
compared to ^^P-MRS^^'^^ approaches, while obtaining full coverage 
of the muscle. As shown in Eq.2, estimation of kf only requires 
measurement of PCr signal at different saturation times. 

Our results are in close agreement with previous ^^P-MR studies. 
Valkovic et al.\ measured a mean /cf in the Gastrocnemius muscle at 
0.31 both at 3.0 T and 7.0 T, using unrealized '^P-MRS. Our 



mean measured A;fin the GL was 0.32 s~^ (Table 1) and 0.31 s~Mnthe 
GM. In an earlier study, Bottomley et al?^ reported a mean /cf of 
0.27 s"^ using a localized MRS method at 1.5 T. Our measurements 
of PCr concentration are well within the range reported in several 
studies reviewed by Kemp et al?^. As seen in Table 1, our measure- 
ments showed significant lower kf in the TA compared to the other 
three muscle groups (i.e. GL, GM, S). These results agree with the 
values we reported at UHF^^, using ^^P-MRI. To the best of our 
knowledge, no ST studies using spectroscopy have been reported 
for that particular muscle. In terms of our measurements, Schar 
et aV reported mean values of 2.27 s in the calf muscle at 3.0 T, 
which is very close to the mean reported values in our work for the 
GL (2.11 s) and the GM (2.15 s). 

The main limitation of our spectrally selective imaging approach, 
compared to ^^P-MRS methods, is the lack of localized measure- 
ments of ATP, which may have a two-fold effect in quantifying /cf 
and Vf. First, in the ST experiments, there could be remaining y-ATP 
signal that is not fully suppressed. If that is the case, there could be 
local errors in the estimation of /cf that are not accounted for by Eq.2. 
However, unlocalized ^^P-MRS measurements in the entire volume 
of the lower leg muscles (Fig. 1) did not show any evidence of unsup- 
pressed y-ATP signal, hence we do not expect that this would affect 
quantification in our study. Second, ^^P-MRS studies assume con- 
stant ATP concentration in the skeletal muscle^ ^ and estimate PCr 
concentration from the ratio of the PCr signal to that of ATP. Our 
study requires the use of external phantoms and in some cases (i.e. 
high BMI subjects) water- content corrections that complicate the 
process of quantifying resting PCr concentration. 

Of note, our imaging method could also be translated for studying 
metabolic activity in the brain. The CK reaction rate has been studied 
with one- dimensional spectroscopy^^ at 3.0 T, while our method 
could provide full-coverage of the brain in similar acquisition time. 

Compared to UHF, ^^P MRI at 3.0 T suffers from a close to three- 
fold decreased signal to noise ratio^^. In our implementation at 
UHF^^, we used a constant saturation time prior to imaging, and 
measured in a separate experiment. However, at 3.0 T, we found 
that mapping was particularly challenging, with the images 
acquired at different inversion times not yielding sufficient SNR in 
order to have accurate fittings of the recovery curve. In the work 
presented here, we used the progressive saturation approach, where 
we saturate y-ATP for different durations. Through a three-para- 
meter fitting of Eq.2 to the data, we recovered Ti and /cf for each 
muscle which was then used in Eq.3 to estimate T/.In this approach. 



Table 1 | Results of the ST experiments in four major leg muscles. Asterisks depict statistically significant differences between muscle groups 
(P<0.05) 

TA GL GM S 



Vol. 



[PCr] 



rr [PCr] 



Vf 



[PCr] 



Vf 



k 



rr [PCr] 



Vf 



(mM) (mM s 



(mM) (mM s 



(mM) (mM s 



(mM) (mM s 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

avg 

SD 



0.29 

0.16 

0.25 

0.28 

0.26 

0.25 

0.23 

0.34 

0.30 

0.27 

0.26* 

0.05 



2.11 

4.12 

2.41 

2.05 

2.42 

2.27 

2.51 

1.69 

1.88 

2.16 

2.39 

0.66 



33.65 
36.60 
29.57 
34.27 
31.12 
33.54 
31.77 
34.22 
29.22 
30.85 
32.48 
2.35 



9.74 
5.91 
7.28 
9.48 
8.03 
8.47 
7.20 
1 1.52 
8.90 
8.29 
8.48* 
1.56 



0.29 
0.26 
0.27 
0.34 
0.27 
0.30 
0.30 
0.37 
0.42 
0.32 
0.32 
0.05 



2.48 
2.29 
2.68 
1.88 
2.74 
2.21 
1.82 
1.50 
1.64 
2.00 
2.11 
0.44 



33.96 
36.31 
32.61 
35.47 
32.91 
34.49 
29.00 
33.22 
28.84 
27.93 
32.48 
2.92 



9.90 

9.50 

8.89 
1 1.95 

8.76 
10.26 

8.84 
12.40 
12.06 

9.56 
10.21 

1.41 



0.35 
0.30 
0.30 
0.29 
0.25 
0.23 
0.30 
0.32 
0.36 
0.37 
0.31 
0.05 



2.00 
2.23 
1.94 
2.26 
2.77 
2.84 
2.1 1 
1.80 
1.81 
1.72 
2.15 
0.39 



36.10 
37.40 
29.96 
34.43 
32.90 
33.76 
30.35 
32.92 
27.83 
28.09 
32.37 
3.25 



12.53 

1 1.33 
9.03 
9.97 
8.18 
7.71 
8.96 

10.65 
9.89 

10.44 
9.87 
1.16 



0.37 
0.30 
0.26 
0.26 
0.29 
0.28 
0.27 
0.30 
0.39 
0.32 
0.30 
0.04 



1.83 
1.99 
2.57 
2.50 
2.45 
2.24 
2.09 
2.03 
1.73 
2.00 
2.14 
0.29 



33.42 
32.81 
27.67 
31.75 
30.37 
31.66 
28.26 
35.17 
26.50 
27.31 
30.49 
2.94 



12.31 
9.75 
7.16 
8.37 
8.82 
9.02 
7.67 
10.71 
10.30 
8.63 
9.27 
1.53 



*significantly lower than GL, GM, S. 
**significantly lower than GL, GM. 
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Figure 4 | Saturation transfer spectroscopy (ST-^^P-MRS) and imaging (ST-^^P-MRI) pulse sequences. A saturation transfer module consists of a 
train of Gaussian pulses, which saturates the y-ATP resonance. Spoiler gradients are used between two consecutive pulses to destroy any remaining 
transverse magnetization. The number of Gaussian pulses defines tsat in each experiment. A) Unlocalized ^^PMRS data were acquired by sampling the free 
induction decay (FID) after implementing the ST module. B) Identical ST preparation module as in ^^P-MRS, followed by a spectrally selective ^^P-MRI 
imaging sequence. 



the PCr signal for all ^sat is between 40-100% of the reference signal 
(Fig. 1), allowing for accurate fitting of Eq.2 to the data. 

In summary, ^^P-MR is unique in its ability to non-invasively 
probe human metabolism in vivo. High-field (3.0 T) clinical magnets 
provide sufficient sensitivity in order to perform imaging experi- 
ments with large tissue coverage and relatively high spatial resolution 
for measuring muscle -specific energy metabolism in the skeletal 
muscle and potential disruptions in local energy metabolism caused 
by disease. 

Methods 

Theory. The CK reaction can be written as: 

PCr-\-ADP-\-H^ XATP-\-Cr (1) 

kr 

with /cf and kj- the pseudo first-order forward and reverse rate constants respectively^. 
One approach for measuring k{ is to saturate the y-ATP resonance for different 
durations (the progressive ST experiment) The saturated magnetization of y-ATP 
is transferred to PCr (Eq.l), resulting in a net decrease of the PCr signal. The forward 
rate k{ can be calculated from the relative decrease of the PCr signal as a function of 
saturation time (tsat)- The CK rate constant A;f multiplied by the PCr concentration 
yields the unidirectional flux of PCr to form ATP. Under fully- relaxed conditions, 
assuming complete saturation of the y-ATP resonance, the magnetization of PCr as a 
function of tsat can be derived from the solution of the Bloch equations, modified for 
chemical exchange (assuming a two-pool exchange system)^'': 

Mitsat) = c[l+kfT,e-^^^'^^'-'] (2) 

Where M(fsat) the magnitude of the PCr signal measured at different times (tsat)> c a 
parameter accounting for direct spill-over effects (i.e. direct saturation of PCr by 
saturating irradiation on y-ATP), and Ti the spin-lattice relaxation time of PCr. By 
measuring the PCr signal for several tsat> can be estimated through a three- 
parameter (i.e. c, kf and Ti) fit of the data to Eq.2. 

The intrinsic spin lattice relaxation (Ti) in the presence of saturating irradiation 
is^^: 

Y7^Y^'^ (3) 



Human Subjects. The study was fully Health Insurance Portability and 
Accountability Act (HIPAA) -compliant and approved by the NYU Institutional 
Review Board. We obtained written informed consent from all participants in this 
study. The methods were carried out in accordance with Food and Drugs 
Administration (FDA) guidelines. We recruited ten non-smoking healthy volunteers 
(seven men, three women, 32.0 ± 3.5 years of age, range 29-41, BMI 22.6 ± 3.9, range 
18.5-30.4) without any medical history of disease affecting muscle function or blood 
flow. We imaged them on a 3.0 T MRI system (Tim Trio, Siemens Medical Solutions, 
Erlangen, Germany) using a dual-tuned (^^P/^H) transmit-receive quadrature 
birdcage knee coil (Rapid MRI, Ohio) (18 cm inner diameter). 

^^P-ST Spectroscopy and Imaging. We performed shimming on the entire volume of 
the lower leg muscles, using the proton channel of the dual-tuned coil and an iterative 
shimming algorithm provided by the manufacturer. We acquired ST spectroscopy 
and imaging data using the pulse sequences shown in Fig. 4. 

Equation 2 holds for complete saturation of y-ATP To confirm the efficiency of 
our ST module to saturate y-ATP, we tested it using an unlocalized pulse acquire 
sequence (Fig. 4A). Acquisition parameters of the sequence are summarized in 
Table 2. We zero-filled data to 8-k data points prior to reconstruction, and applied 
baseline correction for each spectrum separately. 

We acquired ST imaging data using the progressive saturation ^^P-MRI sequence 
shown in Fig. 4B^^. Acquisition parameters of the sequence are summarized in 
Table 2. In both the ST spectroscopy and imaging acquisitions, we used a ST module 
consisting of a train of Gaussian pulses (each 50 ms long, with 56 Hz bandwidth and 
0.85 |iT strength at peak), which saturates the y-ATP resonance. Between two con- 
secutive pulses a 7 ms delay allows for the use of spoiler gradients (magnitude 20 mT 
m"^; duration 5 ms) to destroy any remaining transverse magnetization. The number 
of Gaussian pulses defines tsat in each experiment. 

The imaging sequence was a centric ordered three-dimensional turbo spin echo 
(3D-TSE) with a frequency selective 90° pulse (16 ms duration, 125 Hz bandwidth) 
that excited only the PCr resonance^^. The effective echo was 26 ms, which resulted in 
minimum contamination from ATP signaP*^. A train of 24 non-selective 180° pulses 
was applied after each excitation, using equal area crusher gradients to remove the 
free inductions decays (FIDs) that were produced by imperfections of the refocusing 
pulses. During the ST imaging experiments, the SAR levels never exceeded 5% of the 
maximum allowed under normal operating mode. 

This /c- space sampling scheme resulted in increased blurring in the imaging dir- 
ection where the ETL was sampled. We simulated the effect of signal modulation 
during the ETL acquisition in order to predict the point- spread- function (PSF) along 
that direction, using an average T2 value of PCr of 365 ms based on previous mea- 
surements^^ which were in close agreement with T2 values at 3.0 T reported in the 
literature^'*. 
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Note: FID = free induction decay, NA = not applicable. 



For absolute quantification of PCr in the muscle we imaged a phantom with 
comparable to the in vivo coil loading, using the same ^^P-MRI sequence and para- 
meters (with the ST module turned off), with the exception of repetition time (TR) of 
60 s. We used the phantom to derive the calibration curve of PCr as a function of 
signal intensity. The phantom consisted of three sealed cylindrical tubes containing 
different concentrations (25, 50 and 75 mM) of inorganic phosphate (Pi). 

Water-Fat Imaging. Fat infiltration, especially in high BMI subjects, can reduce 
the volume fraction of lean muscle in a certain volume of tissue affecting 
quantification of PCr concentration. In order to account for the fat content in the 
muscle, we used a three-dimensional gradient echo sequence with three echo-times 
(TE: 2.1, 2.8 and 3.7 ms; flip angle: 3°; TR: 12 ms, total acquisition time: 3 min), with 
the same FOV and frame of reference as the PCr images. The matrix size was 128 X 
128 X 40. We reconstructed water and fat fraction images using the Hierarchical 
IDEAL method^^ on all ten healthy volunteers. These fractions represent the proton 
density fat and water fractions and can be converted to volume fat and water fractions 
as shown previously^^. 

Data Analysis. For each muscle group, we estimated /cf from Eq.2 by fitting the mean 
signal intensity of the segmented PCr datasets at different tsat (excluding slices 1 and 8 
at the edge of the coil). We estimated Vf from the product of /cf and PCr concentration 
in the same volume of the muscle. We compared /cf and Vf values among four different 
muscle groups [Gastrocnemius Lateralis (GL) and Medialis (GM), Soleus (S), and 
Tibialis Anterior (TA)], with paired t-tests with a 5% significance level. 
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